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Abstract. Literature recorded that articles that used weighted goal programming in pre-1990 application 
papers is 25% whereas its use in the period 1990-2000 is 41%. The reasons for this increment are due to the 
greater flexibility by the weighted goal programming, and the aim of the decision makers to do more ‘trade- 
off analysis and direct comparison between goals. An efficient solution algorithm for solving weighted goal 
programming problems is presented. The al gorithm incorporates rigid constraints in the model and utilized 
variant of simplex method. A step by step illustration of the procedure is provided. The computational 
method is efficient. 


1. Introduction 

In 1948, Dantzig [1] developed a systematic procedure for solving linear 
programming problems (LP) called simplex method which is efficient in practice and is a 
great improvement over earlier methods of solving LP such as Fourier-Motzkin elimination, 
even though it is exponential in complexity. Lp had received substantial international 
exposure and attention, and was hailed as one of the major developments of applied 
mathematics. Today, linear programming is the most widely employed of the methods used 
by those in such field as operations research and management sciences. However, as with 
many quantitative approaches to the modeling and solution of real problems, LP has its 
blemish, drawback and limitations. Its inability to directly and efficiently handle problems 
involving multiple objectives and goals: thus restricting the users of the technique to 
narrowing their problems to a single objective function. But in real life, decision situations 
frequently are characterized by multiple and conflicting objectives. For instance, Rifai [2] 
explained that, besides profit maximization, management has many other objectives to 
attain, including improved market share, maintaining full employment, providing quality 
ecological management, improved relationships with the public, production of a quality 
product, minimized pollution and the realization of a certain rate of return in investment. 
Many of these objectives may conflict, and so compete for scarce resources. LP is an 
inadequate technique for solving problems of this type. One way of handling problems with 
multiple objectives is to choose one of the goals as the supreme goal and to treat the others 
as constraints to ensure that some minimal ‘satisficing’ level of the other goals is achieved. 


2010 Mathematics Subject Classification: 90C29. 
Key words and phrases: weighted goal programming problem. 
©2013 Science Asia 
1/11 


T UKAMAKA CYNTHIA ORUMIE & DANIEL EBONG 

As a result of this drawback, a number of programming methods have evolved based on the 
Dantzig algorithm to solve programming problem of this situations. One of such class of 
problems is referred to as goal programming problems. 

Goal programming (GP) can be considered as a popular branch of multi-objective 
optimization which itself is a part of multi-criteria decision analysis. This can also be 
considered as a generalization or extension of linear programming, in which more than one 
and in general conflicting objectives are to be optimized. GP technique was initially 
developed by Charnes and Cooper [3], Lee [4], Ignizio [5] and several others to handle 
multiple objective situations that cannot be addressed by linear programming. The essence 
of the technique is the achievement of the “best possible” solution, which comes as close as 
possible to meeting the goals. 

However, Olson [6], described lack of an algorithm capable of reaching a 
compromised solution in a reasonable time as a major setback in applying GP.(see also 
Schniederjans and Kwaks [7] , Baykasoglu [8]. The most commonly used goal programming 
solution methods were introduced by Lee and Ignizio based on Dantzig's simplex method 
(see Authur and Ravindran [9], and Schniederjans and kwaks. Both methods require 
columns in the simplex tableau for positive and negative deviational variables. They also 
require separate objective function rows for each priority level, all of which add immensely 
to the computational time of the solution method. Tamiz and Jones [10] pointed out that a 
faster mathematical programming based solution method is given by Ignizio [11]. This has 
the same effect as imposing constraints on the unwanted deviational variables but reduces 
the problem in terms of number of variables and size of feasible region while keeping the 
number of constraints constant. Orumie and Ebong [12] developed an efficient method of 
solving a generalized linear goal programming problem. 

However, weighted goal programming and Lexicographic goal programming are the 

most popular variants of the GP model as shown in Alp et al [13] and many others. The two 
methods do not generate the same solution and neither is one method superior to the other 
because, variantis designed to satisfy certain decision makers' preferences. 
The purpose of this research is to present an efficient solution method for strictly linear 
weighted goal programming problems. In section two, the linear weighted goal 
programming problems will be formulated whereas the statement of problems and 
methodology will be presented in section three and four respectively. The new algorithm 
for solving weighted goal programming will be formulated in section five, followed by step- 
by-step illustration of a case problem in section six. Section seven and eight presents the 
conclusion and further research respectively. 


2 Weighted Goal Programming (Non-Pre-emptive) 


WEIGHTED GOAL PROGRAMMING PROBLEM — 
If the decision maker is more interested in direct comparisons of the objectives, 
then weighted or non pre-emptive goal programming (WGP) should be used. This variant 
allows for direct trade-offs between all the deviational variables to be penalized by placing 
them in a weighted, normalised single achievement function. I.e all the unwanted 
deviations are multiplied by weights, reflecting their relative importance, and added 
together as a single sum to form the achievement function.( Weights are attached to each of 
the objectives to measure the relative importance of deviations from their target). I.e It is 
assumed that all the goals are of equal important with only slight differences that can be 
measured by assigning weights to these goals. WGP handles several objectives 
simultaneously by establishing a specific numeric goal for each of the objectives and then 
find a solution that comes close to each of these goals. The overall purpose is to minimize 
the weighted sum of deviations of these objective functions from their respective target 
values. It is important to realize the fact that deviational variables measured in different 
units cannot be added because of the principles of incommensurability. Hence each 
unwanted deviation is multiplied by a normalization constant to allow direct comparison. 
Popular choices for normalization constants are the goal target value of the corresponding 
objective (hence turning all deviations into percentages) or the range of the corresponding 
objective (between the best and the worst possible values, hence mapping all deviations 
onto a zero-one range). 
The algebraic formulation of a WGP as presented in Min and Storbeck [14] is given as: 


m 


min z= (w, d, +w,d,) (1) 

s.t 

Zax, +4, -d; =b, (i =1,2, m), (2) 
3J 

S58, dy 20,w,>0 (3) 
(i =1,2,..,m: jf =1,2,3...,n) (4) 


where w, ,w, are non-negative constants assigned to negative and positive deviational 


variable d, and d, (= 0) respectively. Equation (1) represents the achievement function 


which measures the unwanted weighted sum of all the goal deviational variables to be 
penalized. Equation (2) is the objective goal whose deviation from the target value bi is to 
be minimized in the achievement function. Tamiz et al.[15] reported that articles that used 
weighed goal programming in pre-1990 application papers is 25%. Tamiz and Jones [16] 
reported that its use in the period 1990-2000 41%. The reasons for this increment are due 
to the greater flexibility by the weighted goal programming, and the aim of the decision 
makers to do more ‘trade-off analysis and direct comparison between goals. 
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3. Statement of problems 
According to Min and Storbeck[14], formulation of a WGP is given in equations (1)- 
(4). The above statement (model) omitted the fact that a goal programming model may 


include a rigid constraints like Ð a,x, < b,.In order to incorporate the rigid constraints in 


the model, the researcher decided to augument the above rigid constraints to equations 
(1)-(4) and to develop new algorithm for linear weighted goal programming as will be 
shown in section five. 


4 Methodology /Initialization (Initial Table formulation) 

The new goal programming algorithm is formulated into initial tableau in the same 
format with that of Orumie and Ebong . The different is that in the new algorithm, the 
deviational variables that appeared in the achievement functions with their weight 
attached to them, together with slacks from the rigid constraints if exists forms the basis 
and the column of the slacks variables from the rigid constraints added to the non basic. 
The table 4.1 below represents the initial table. The decision variables xj, the slack variables 
from the rigid constraints, and the deviational variable that appeared in the achievement 
function are placed in row 1, column2. Last Column constitute the right-hand-values, bi. 
The coefficient of the decision variables aij, the slack variables (s) from the rigid constraints, 
and the deviational variable that appeared in the achievement function (ci) which 
constitute the coefficient matrix G =[(¢, )], where gy= f{aj, s, cyM} > 


8, =a aa 6 AnS CaCa » - €,+) be the i row of G and g the jè column of G are 

placed in column 2. The weihted factors wefor the deviational variable, including slacks 

variables are presented in column 1 which forms the basis. 

It starts by not including the deviational variables column that did not appear in the basis 

from the tableau ( ie while searching for the optimal solution), but others can be 

augmented when necessary since positive deviational variables columns coefficient is the 

same as negative of the negative deviational column coefficient as represented in table (4.1) 
below. This is because d= -d{+. The new algorithm starts with the highest weighted row 
constraints. 

Let di” = ditor dy €z> v= F E Oe 


|- if v is positive sign 


) 


Let c0 be the coefficient of dy € z in the basis, then initially, the coefficient of a‘ is given 


by 
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gyei 1 if v is negative sign 
l |-1 if v is negative sign 


Here a” 


are the deviational variables in the objective function, and t is the number of 
deviational variables in the objective function. Let g, be the coefficient of weighted row(s) 


corresponding to wxwith initial table below. 


(Table 4.1) Innitial Table of the new algorithm 
Variable in basis with wx.Cp X X2 m XnS dim dom. . deo Solution value b; (rhs) 


. Ain $1 C110 C110). . . CHM 


. aan S2 C21 C220. ..Car™ 





5 The New Algorithm For Weighted Goal Programming (WGP) 


Consider the weighted goal programming where weights are assigned to the goal 
constraints, (deviational variables in the achievement function) .The WGPP formulation for 
n variables, m goal constraints, t deviational variables and F weight factors that that 
includes the rigid constraints is defined below: 


r 
min z=} (w,d, +w,d}) fr fci, c{1,2,. . m} (5) 
f 


such that 


ax +d -d`=b 
> y ri l 


i i 


(6) 
ax, <b 
l (7) 
Where 
w, 20, md, a, 20 x,,d;,d} 20 (8) 


for CHllsnm’ j=1,2,3.. n) (3) 
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where w= f” weighted factor f= 1,2,...F 
(wid, + wd; ) areset of the weighted sum of the deviational variables in z. 


e.g 
Min (4a2++3d) + 2(d3+d#)) 
st x 43x. <120 
2X1 +x2 + 2x3 <80 
10x%14+15x24+d -di+= 100 
Xi + 2x2 +dz-d+ = 50 
2x1 + x + x4+d3 -a3+ = 100 
Xi tdy-dt+ = 80 
Xi, X, X, X4, dv, dit, dz, dot, dz, dst , dy „dat > 0 

f=1,2,3, j=1,...4, i=1,2, ..6, t=1,2,3,4 


The algorithm (WGP): 


Let w) =w, or w? bethe f” weight assigned to a’ ,then, 


Step 1. Initialization: 
Seth & 1 
Step 2. Feasibility: 
If b, = 0 foralli=1,2,. . . m, go to Step 8. {solution optimal} 
Set b, < le,| foral I =1,2,. . .,m, {ensure feasibility} 
Step 3. Optimality test: 


If g, <0 forallj+ pivot colum ,hei, go to Step 7. 


{all coefficients of weighted row, h non positive ,so w} :v = + or — is satisfied} 


Step 4. Entering variable: 
Entering variable is the variable with highest positive coefficient in the row 


g, hei, for the w} (d; +d; ) rows of the objective functions which does 


not violate the satisfied weighted deviational variable. If g,,h  {1,2,...F} is 


the highest coefficient, but has been previously satisfied or has more weight 
attached to it than the leaving variable under consideration then consider the 
next higher value on the same row, otherwise go to step 7. 

(The weight attached to the entering variable should be placed alongside 
with it into the basis). 


Step 5. 


Step 6. 


Step 7. 


Step 8. 
i) 
ii.) 
iii.) 


Note: 
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In case of ties {g,,....g,,,}, then the entering variable is the variable for 


. ae ‘Tis; : 
which min 3“: :g, >0$ ismaximum. j,:p =1.... 0 
h [ E Jp J 


{ties in the weighted rows}. 


Leaving variable: 
If y, is the column corresponding to the entering variable in step 4, then the 


leaving variable is the basic variable with minimum 
b . 
We g > OSI n . {e , is the pivot column}. 
g a c20 “Yo 


In case of ties, the variable with the smallest right hand side leaves the basis. 
Interchange basic variable with non basic: 
Perform Gauss Jordan row operations to update the table.If w” is still in the 


basis (Cp), go to Step 3. 

Increment the process: 

Sete «— f +1 lf < F, go to Step 3. Satisfied weighted variable will not 
reenter for the lesser one to leave, instead the next highest coefficient 
variable enters. 

Solution is optimal: 

The coefficient of the weighted rows are all negative or zero 

The right hand sides of the weighted rows are all zero 

The weighted rule is satisfied. 

The optimal solution is the value of wi (d, +d;) in the objective function as 


appeared in the last iteration. 
Just as in the method of artificial variables, a variable of higher or equal 
weight that has been satisfied should not be allowed to re-enter the table. In 


this case the next higher coefficient of ¢,, will be considered. 


6 Step by step illustration of the algorithm on weighted goal programming problems. (The 
algorithm WGP) 


(i) min z= (2d, +d2+) 


S.t 4x +8x, 


+d - d? =45 


8x1 + 24x2 + d} -d} =100 


X1 + 2x25 10 
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Source: Cohon, [17] 


[4 8 1 0 0 0] 
The G matrix of the above problem is | $ i ! ° J with the bi 
l1 2 0 0 1 0 | 
li 0 0 0 0 | 
[45 ] 
_ | 100 | 
lio | 
o |] 


Table (6.1) below is the initial table of the problem (i). The first column is the column of 
the deviational variables to be minimized with weight attached to them and the slacks from 
the hard constraints (basic variables). sı and s2 are the slacks variables attached to the 
third and forth constraints respectively (dummy variables) . 


Table (6.1) Initial table of problem i 





x2 enter 
d, leaves 
Step 1. Set f =1. 
Step 2. Jb, #0 for i=1.,.., m £ feasible.} 
Step 3. 3 g,,>0 fr for some j . {Solution not optimal} 
Step4. Max{g,,} = max{4,8,1,0,0,0}= 8 at ¢,, .So,x2 enters. 
Step5. Min (ee: ga > 0}= min (4,25) = a + . dz* leaves the basis. 
Step6. Perform the normal Gauss Jordan row operations to update the new table 


(see table 6.2) and checkif w, is still in the basis (Cs) to test for optimality. 


Table (6.2) 1st Iteration for problem i 
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e el see 


25/6 


Si leaves 





In table (3.5.2), w, is not 


satisfied, so we go back to step 3. 


Step3. 3g,,>0 fr for some j . {Solution not optimal} 

Step4. dmax{ g,,} = Max{+,0,1,1/3,0,0} = 4/3 at g,, X1 enters. 

Step5. Min = ga> 0}= min {#, 4,5,9} =5 at A .S1 leaves. 

Step6. Perform the normal gauss Jordan row operations to update the tableau (see 


table 6.3) and check if w, is still in the basis (Cs) to test for optimality. 


Table (6.3) 2"4 Iteration for problem i 





In this table, w, is not satisfied, so we go back to step 3. 


Step3. g,, <0 fr all j . {Solution optimal} 

Step4. dmax{ g,,} = Max{0,0,1,0,-4,0} =1 at dz* but recall that higher weighted 
factor will not re-enter the basis for lower one to leave. So analyses stop. 

Step8. Thus, table (6.3) is therefore an optimal tableau with z=10,x1=5,x2=2.5,d1- 
=5 and s2=4. 


7 Conclusion 
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The new method is applied to various problems of different variables sizes, goals 


and constraints. The proposed method is efficient in reaching optimality and its 
formulation represents a better model than the commonly used. 


8 Area for Further Research 
Further research should be carried out on the modification of the new algorithm to 
solve other types of GPP. 
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